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ABSTRACT 

We present a combined analysis of the low-mass initial mass function (IMF) 
for seven star-forming regions. We first demonstrate that the ratios of stars to 
brown dwarfs are consistent with a single underlying IMF. By assuming that the 
underlying IMF is the same for all seven clusters and by combining the ratio 
of stars to brown dwarfs from each cluster we constrain the shape of the brown 
dwarf IMF and find it to be consistent with a lognormal IMF . This provides the 
strongest constraint yet that the substellar IMF turns over oc M~°, a < 0). 

Subject headings: stars: initial mass function — pre-main sequence — forma- 
tion; brown dwarfs 



1. Introduction 

Speculations concerning the exist ence and frequency of brown dwarfs ca n be traced to 
before the introduction of the term (IKumaii Il963l : iHayashi fc Nakand 119631 ). Since then, 
wide-field surveys have uncovered h undreds of candida tes in the field and revealed two new 
spectral types, the L and T dwarfs (IKirkpatrickl 120051 ) . Yet the frequency of brown dwarfs 
compared t o star s has remained a topic of confusion and debate. In a pioneering work, 
Reid et al.l (119991 ) attempted the first census of the substellar initial mas s function (IMF) 
based on results from the Two Micron All Sky Survey (ISkrutskie et al.ll2006r). They pr esented 



evidence for a low-mass IMF that was more shallow than a Salpeter (jSalpeter 



suggesting that brown dwarfs were not a significant contributor to dark matter. 



9551) slope. 



Allen et al. 
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( 120051 ) used a Bayesian approach to constrain the power-law slope below 0.08 Mq to be in 
the range —0.6 < a < 0.6 with a confidence level of 60%, where a Salpeter slope is a = 2.35. 
These results indicate that, although brown dwarfs do not contribute sig nificantly to th e 
mass of typical stellar populations, they might still be as abundant as stars (IChabrierll2002l ). 



The classical approach to deriving the mass function for stars and substellar objects is 
to take an observed luminosity function and apply a mass-luminosity relationship in order 
to derive the present-day mass function. Then, corrections, based on the theory of stellar 
evolution, perrn it one t o estimate an ini ti al mass fu n ction from the present-day mass function 



(see e.g. IScald (Il986l ): iKroupal (1200 ll ): IChabrierl (120031 ) for complete descriptions of this 
process). The confounding variable in these analyses is the star formation history of the 
Galactic disk, which is vital for substellar objects whose mass-luminosity relationship evolves 
with time. 

A different approach is to use star clusters of known age as laboratories to measure 
the IMF. Open clusters are in principle good candidates because of their richness. Yet 
they suffer from th e effects of dynamical evolution, mass segregation, and evaporation (e.g. 
Lada fc Ladal 120031 ). Young (< 10 Myr) embedded clusters are attractive alternatives as 
they are compact and rich (from hundreds to thousands of stars within 0.3-1 pc), and yet 
to emerge as unbound OB/T associations, and the low mass objects are 10-1000 times more 
luminous than their older open cluster counterparts (0.1-16 Gyr) because they shrink and 
cool as they age. 

Indeed, embedded clusters have been the targets of aggressive photometric and spec- 
troscopic su rveys in an a t temp t to search for variations in the IMF as a function of initial 
conditions. iMeyer et al.l (120001 ) found that the ratio of high- mass (1-10 Mq) to low-mass 
(0.1-1 Mq) stars for an ensemble of young clusters within 1 kpc was consistent with (1) each 
other and (2) having been drawn from th e field star IMF. Mor e recent studies have pushed 
well into the substellar mass regime (see iLuhman et al.l (120071 ) for a recent review) . There 
have bee n some claims for va riations in the brown dwarf IMF between nearby star-forming 
regions. iBricefio et al.l (120021 ) argued that the low-density Taurus dark cloud had a dearth 



of brown dwarfs compared to the rich Orion Nebula Cluster (ONC). However, this prelimi- 
nary result has been updat ed as additional data have become avai lable and as the statistics 
improved for both clusters (iGuieu et al.ll2006l : ISlesnick et al.l 120041 ). 



Here we use observations of seven nearby star clusters to constrain the combined brown 
dwarf IMF. In section 2, we describe the data, illustrate that there is no strong evidence for 
variation in the substellar IMF between the star-forming regions, and outline our approach 
to constrain the low-mass IMF. In section 3 we present our results, and in section 4 we 
discuss our results in the context of previous work as well as theories of star (and substellar 
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object) formation. 



2. The Approach 



We have compiled the ratio of stars to brown dwarfs in nearby, weU-studied young 
embedded clusters and the Pleiades. The regions included in this study are described briefly 
below, where the ratio of stars (0.08-1.0 Mq) to brown dwarfs (0.03-0.08 Mq) is calculated. 
For all the regions, we consider the system IMF, uncorrected for multiplicity within 200 
AU. The sample is focused on embedded clusters, in which spectroscopy has been used to 
determine the age of the cluster, field star contamination has been taken into account, and 
an extinction-limited sample has been defined. Furthermore, we have included the Pleiades, 
because it is one of the best-studied open clusters and bacause its substellar IMF has been 
estimated. The break poin t at 0.08 Mq has been adopted in accordance with the break 
point for the iKroupal (120011 ) IMF, similar to the characteristic mass in the IChabrierl (120031 ) 
single object IMF. Only a few of the clusters adopted here have the IMF derived in an 
extinction-limited sample reaching 0.02 Mq and we have opted for 0.03 Mq as a lower mass 
limit to obtain a larger sample of clusters. 



Taurus. iLuhmanI (120041 ) imaged a 4 deg^ region of Taurus that focused on the denser 
filaments, to identify cluster candidates. Candidates were confirmed as cluster members,by 
use of follow-up intermediate-resolutionoptical spectroscopy, on the basis of their effective 
temperature, luminosity, and spectral features. In total, 112 objects were confirmed members 
with derived masses between 0.03 and 1.0 Mq and extinctions Ay < 4 mag. Some 96 objects 
were stars and 16 were brown dwarfs. Thus, the ratio of stars to brown dwarfs in Taurus 
was found to be i? = 96/16 = 6.0^*^2 where the errors are estimated using the method of 



Gehrels I (119861 ). 



IC 348. iLuhman et all (120031 ) imaged a 42'x28' region of the IC 348 cluster to iden- 
tify cluster candidates. By the use of intermediate-resolution spectroscopy, most of the 
candidates were confirmed as cluster members, on the basis of their effective temperature, 
luminosity, an d spec tral features, which indicated that the objects were young. In total, 
Luhman et al.l (120031 ) found 168 cluster members with masses between 0.03 and 1.0 Mq and 
extinctions Ay < 4 mag. The ratio of stars to brown dwarfs was found to be i? = 8.3] 



1-3.3 
-2.6- 



Mon R2. lAndersen et al.l (120061 ) imaged the central I'xl' of the embedded cluster asso- 
ciated with Mon R2 by utilizing the Near-Infrared Camera andd Multi-Object Spectrometer 
on board the Hubble Space Telescope [HST). An extinction-limited sample Ay < 10 mag was 
defined and a total of 19 objects were detected with masses between 0.03 and 1 Mq. The 
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ratio of stars to brown dwarfs was found to be i? = 8.5 



+13.6 
-5.8 • 



Chameleon 1. iLuhmanI (120071 ) obtained an extinction- limited sample in Chameleon 1 
that was complete down to 0.01 Mq for Ay < 5 mag, by use of observations of a 0.22° xO. 28° 
region with the Advanced Camera for Surveyes on board HST and a subsequent spectroscopic 
follow-up of cluster member candidates. The sub-sample from 0.03 to 1 Mq includes 24 
objects and the ratio R was found to be i? = 4.0^^'^ 
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Pleiades. The Pleiades is one of the best-studied open clusters, and numerous deriva- 
tions of the IMF have been published. Here we focus on the survey by iMoraux et al.l (120031 ) 
who covered a 6.4 deg^ region of the Pleiades. The survey had a saturation limit of 0.48 
Mf^. For higher rn a sses, the survey was combined with a mass function built using the 
Prosser &: Stauffer I (Il998l ) database. The Pleiades suffer relatively low {Ay < 1 mag). 



mostly uniform, extinction, with negligible impact on the completeness of this sample, so 
we did not apply a reddening criterion. The ratio of stars to brown dwarfs was found to be 
R = 4:.9tl:l 



The Orion Nebular Cluster. The ONC has been the subject of extensive stu dies ( jHillenbrand 



19971 : iHillenbrand fc Carpenterll2000l : iLuhman et al.ll2000l : iMuench et al.ll2002r) . We ta ke the 



adopted ratio of stars to substellar objects from the study of ISlesnick et al.l (120041 ). The 
total sample, covering the central 5.1'x5.1', contains approximately 200 objects with masses 
between 0.02 and 0.6 Mq and Ay < 15 mag. Using their Figure 14, and extrapolating the 
slope from 0.08-0.6 to 1.0 Mq (one additional bin in their plot), we arrive at a ratio of stars 
to substellar objects oi R = 3.3^0;®. 



NGC 2024. The ratio of stars to brown dwarfs in NGC 2024 was found by lLevine et al. 



(120061 ) from their photometric and spectroscopic study, covering the central lO'xlO'. They 
assigned masses to the photometric objects on the basis of th e mass distribution in each 
magnitude bin, determined from the spectroscopic sample as in ISlesnick et al.l (120041 ). The 
result was that a total of 148 objects in their survey area has masses between 0.02 and 1 Mq 
and extinctions Ay < 15 mag. Based on their Figure 9, we find that there are 27 objects 
between 0.03 and 0.08 Mq resulting in a ratio of stars to substellar objects of R 



3.8l?i. 



Table [T] shows the ratio of stars to brown dwarfs for nearby embedded clusters and the 
Pleiades, as described above, and the distribution of ratios is shown in Fig, [H The weighted 
mean of the ratios is found to be 4.3, and the standard deviation of the weighted mean is 
1.6. All of the measurements presented are consistent with the weighted mean within 2a. 
There is thus little evidence for variation in the low-mass IMF between the different regions 
and we have adopted the hypothesis that the IMF is universal. Under this assumption, the 
complete set of IMF determinations can be combined to place constraints that are stronger 
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than for each of the individual measurements. 



The Results 



For each cluster, we have calculated the probability of obtaining the observed ratio of 
stars to brown dwarfs for a given IMF or greater. The ratio of stars to brown dwarfs drawn 
from a given sample size with an assumed IMF is determined by the bi nomial th e orem . 



The predicted distribution of ratios from both segmented power-laws and a IChabrierl (12005 

dN ^ ^ (log m-log mo)2 
dlogm ^ 2a^ 



mo 



0.25, cr = 0.55) lognormal IMF for a cluster of 100 objects 
with unresolved binaries is shown in the lower panel in Fig.[Tl The peak mass in the lognorma l 
is slightly higher, and the width is slightly more nar r ow th an is presented in IChabrierl (120031 ). 
The chan ge in the best-fit parameters in IChabrieii (120051 ) is due to an updated luminosity 
function (IReid et al.l |2002| ). A similar increase in the peak mass has been suggested by 
(ICovev et allboosh . 



The slope of the segmented power-law between 0.08 and 1.0 Mq was chosen to be 
a = 1.3, and the slope has been varied b elow 0.08 Mg^ in t he range —0.6 < a < 0.6, which is 
the 60% confidence interval presented by I Allen et al.l (120051 ). It is clear that the rising and fiat 
IMFs (a = 0.6, and 0.0, respectively) are difficult to reconcile with the observed distribution 
of ratios. We have quantitatively assessed the likelihood of obtaining the observed ratios 
from an assumed IMF as follows. For each of the seven measurements, the probability of 
obtaining that ratio or higher, assuming an underlying IMF, is calculated by adopting the 
binomial theorem. The product of the seven probabilities is then calculated. We find these 
values, which we refer to as the binomial tail product, or BTP, to be 0.0012, 2.2 x 10~^, 
1.8 X 10~^^, and 1.0 x 10~^^, for a Chabrier, falling, fiat, and rising IMF, respectively. If 
each cluster sample was drawn from the assumed underlying IMF, and if each cluster had 
an infinite number of objects, we would expect the combined product of this statistic for a 
sample of seven clusters to be 0.5'' = 7.8 x 10~^. The lognormal IMF appears to reproduce 
the observed ratios best, followed by the falling power-law IMF. 

How consistent are the measured ratios with a Chabrier IMF and with what confidence 
can other IMFs be ruled out? We have investigated that question by performing Monte 
Carlo simulations. We created an artificial set of seven clusters, each containing 100 objects 
(the median number of objects in our sample). The 100 objects are then assigned masses 
according to the assumed underlying IMF, and the ratio of stars to brown dwarfs for each 
cluster is determined. For each of the ratios, the probability of observing that value or higher 
is calculated and the seven probabilities are multiplied, as was done for the observed set of 
clusters. The BTP for the observed clusters is then compared with the distribution of BTPs 
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just derived. Because each factor in the BTP is drawn from a binomial distribution (of 
varying shapes), each IMF gives the same expected distribution of BTPs. Figure [2] shows 
the cumulative distribution of BTPs for a set of 10,000 simulations. 

Overplotted are the probabilities obtained above for the observed set of clusters as- 
suming the four different underlying IMFs. We find that 37% of the simulations have a 
probability equal to or lower than what was found assuming a Chabrier IMF, and in only 
~0.05%-0.1% of the simulations is the probability equal to or lower than found assuming a 
falling power-law IMF. In none of the simulations did the low probabilities for the flat or 
rising power-law IMFs occur (P < 0.01 %). The results indicate that the IMF is falling in 
the brown dwarf regime and that the Chabrier IMF is consistent with the observations. 



Discussion 



The results on the IMF presented here are based on the system IMF, including binaries 
unresolved within 200 AU. As such, they may b e difficult to compa re directly with the locally 
derived (within 20 pc) field IMF discussed in lAUen et al.l (120051 ) that suffers from a much 
smaller fraction of unresolved binaries. Yet the overall binary f requency for ultra-cool dwarfs 
(M6 and later) appears to be low (~20%. lBurgasser et al.ll2007l ). and furthermore the relative 
number of companions with separ ations >15 A U and mass ratios q > 0.4 may be extremely 
low around very cool stars ~ 1%; lAUeru (120071 ). 



Indeed, if the companion mass ratio distribution follows the Chabrier IMF at wide 
separations, then one could expect f ewer very low mas s companions as one surveys pro- 
gr essively lower mass primaries (e .g. ISiegler et al.ll2005l ). consistent with the observations 
by [McCarthy &: ZuckermanI (120041 ) . If the IMF follows a Chabrier IMF in the brown dwarf 
regir ne below 0.03 MfD (say, down to the opacity limit for fragmentation of ~ 0.001-0.004 
Mq] IWhitworth fc Stamatellod (120061 ) . then the number of stars below 1 Mq will outnumber 
brown dwarfs 4.7 to 1. 

The sense of our results, that the mass function is falling in th e BP regime, i s consi stent 
with various ideas put forward to explain th e shape of the IMF (IBonnell et al.l (120071) and 
references therein). Building on the ideas of Larson ( 2005 ). Bonnell et al. ~ (j2'o06 ) produced 
an IMF that is only wea kly dependen t on the Jeans mass through dynamical interactions 
in the cluster. Howev er, lAlleru (120071 ) show that the turbulent fragmentation rnodels by 
Bate fc Bonnelll (120051 ) predict too few low-mass binary systems. iGoodwin et al.l (120041 ). on 
the other hand, suggest that the IMF should peak at higher masses in regions with low 
turbulence (e.g. Taurus) which would result in a higher ratio of stars to brown dwarfs. The 
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lack of a strong variations in the ratio of stars to brown dwarfs is a problem for the turbulence 
models in general; for example, magnetic turbulence mod els predict strong variatiori s in the 



low-mass IMF as a function of Mach number and density (jPadoan &: Nordlundll2002l ). If the 



preliminary results indicated here are borne out through further observa tions, then models 



that d e pend only weakly on initial c onditions would be required (e.g. lAdams fc Fatuzzo 



fll996f ): lHennebelle fc Chabriei] tOOm . 



Possible IMF variations at least within 1 kpc are smaller than can be detected by 
comparing the currently observed clusters. Thus, there are two challenges in detecting IMF 
variations: (1) One needs clusters with a well-sampled population to minimize the inherently 
stochastic nature of populating an IMF, and (2) a larger set of clusters is needed to detect 
even small IMF variations with initial conditions. Although it appears that the variations 
in the IMF down to 30 Mjup are modest, we still expect that variation s will be seen at the 



lowest masses where the opacity limit for fragmentation can be reached (ILow fc Lynden-Bell 



19761 ) and the metallicity of the star forming region could be imprinted in the lower mass 



limit. 
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Table 1. Ratio of Stars to substellar Objects in Young Clusters. The distance, age, number of objects in the sample, 
and the extinction limit used for the embedded clusters are given. The four last columns gives the probability of the 

observed ratio having been drawn from the assumed IMFs. 
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Fig. 1. — Top panel: Histogram of the observed ratios of stars to brown dwarfs described 
in the text and summarized in Table [TJ Bottom panel: Binomial distribution for a cluster 
with 100 objects drawn from either the Chabrier [solid line), the falling [a = —0.6, dotted 
line), the flat [a = 0; long-dashed), or the rising [a = 0.6; long-dash-dotted line) IMF. 
Distributions that continue to rise in linear mass units below the hydrogen burning limit are 
least consistent with the observations. 
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Fig. 2. — Test of the distribution of the product of probabihties if seven clusters are randomly 
drawn or higher from a Chabrier IMF. For each the probability of obtaining the observed 
ratio of stars to brown dwarfs is calculated and the product of the seven probabilities is 
determined for each of the 10 000 simulations. The vertical lines indicate the combined 
probability of obtaining the observed ratios of stars to brown dwarfs for the Chabrier IMF 
{right vertical dotted line) and the power-law IMF that is falling in linear units in the brown 
dwarf regime {a = 0.6; left vertical dotted line). The probabihties for the flat and rising IMF 
are both outside the plotted range and did not happen in any of the Monte Carlo simulations. 



